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Abstract 



Phase separation dynamics with an initially non-uniform concentration 
are studied. Critical and off-critical behavior is observed simultaneously. A 
mechanism for an expanding phase separated region is demonstrated and the 
time dependence of the concentration is determined. The final equilibrium 
state consists of a planar interface separating one phase from the other. The 
evolution to this state is characterized by an experimentally observable flux, 
j, crossing this interface. We find that j ~ i~ 2 / 3 if patterns are formed in 
the bulk and j ~ i~C 2 if the bulk remains homogeneous. The results are 
explained in terms of scaling arguments which are confirmed numerically. 
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Phase separation dynamics of a binary mixture is a prototype nonequilibrium process 
exhibiting features such as domain growth, front propagation, pattern formation and pattern 
selection [1|-|3| . In a typical experiment, a binary mixture in the homogeneous phase is rapidly 
quenched to inside the coexistence curve. For a near critical quench convoluted domains 
of each component are formed, while in an off-critical quench, droplets of the minority 
phase form within the majority phase background. The phase separation process exhibits 
dynamical scaling and the characteristic domain size, R, usually grows as a power law in 
time. If the dynamics conserves the order parameter and hydrodynamics is irrelevant, R(t) 
grows as t 1 / 3 at late times ||,|5|]. 

Usually critical and off-critical quenches occur separately. However, if the initial state 
is inhomogeneous, these processes can be concurrent. For example, in the binary liquid 
experiments of Jayalakshmi et al., the fluid was prepared with a linear gradient concentration 
profile ||. A spinodal decomposition region was observed in the middle of the cell with 
slowly expanding nucleation regions on either side of the central region J7|. Kolb et al. 
performed simulations of the kinetic Ising model using both a linear gradient profile and a 
step function profile as the initial states ||. Their results were interpreted using gradient 
percolation theory and they concluded that the concentration profile does not change with 
time. 

In this letter we study phase separation with inhomogeneous initial conditions via a 
numerical integration of the time-dependent Ginzburg-Landau equation for a conserved field 
(the Cahn-Hilliard equation) ||. We discuss an additional mechanism for the expansion, 
with time, of the area in which phase separation has occurred. We show that the kinetic Ising 
model simulations || describe early time behavior and follow the change in the concentration 
profile with time. However, since our method is very effective in obtaining the long time 
behavior ||, our emphasis, will be on the late stage dynamics not discussed previously. To 
reach the final equilibrium configuration of the two co-existing phases separated by a planar 
interface, mass must be transferred from one side of the system to the other. We find that 
the mass current density, j, behaves as j ~ t -2 / 3 if patterns are formed in the bulk after the 
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quench and j ~ t -1 / 2 if the bulk remains homogeneous. The results are explained in terms 
of scaling arguments based on, in the first case, the formation of droplets, and, in the latter 
case, diffusion. We confirm the scaling assumptions numerically. 

We describe the phase separation dynamics using the Cahn-Hilliard equation, 

where c(r, t) is the deviation of the concentration from its critical value, F[c] is the Ginzburg- 
Landau c 4 free energy. The local chemical potential, /z, is then 

Mr, t) = = -c(r, t) + c(r, t) 3 - V 2 c(r, t). (2) 



Here c, r and t are scaled to set all coefficients to unity fllPR . We assume a deep quench so 
that an additive noise term in Eq. ([!]) can be discarded. Therefore, any phase separation 
must be triggered by randomness in the initial condition ||. 

We numerically integrated Eq. (|l|) using an Euler discretization on an L x x L y lattice 
with mesh size Ax = Ay = 1.0 and time step At = 0.01. Periodic boundary conditions are 



used in the x-direction and reflecting b.c. in the y-direction ('cylindrical symmetry') |p"T 
Two types of initial conditions are examined. We consider a linear gradient profile, 

c(r,0)=c -^/ + £(r), (3) 

Ly 

where r = (x,y), cq is the absolute maximum of c and £(r) is an uncorrelated random 
variable with standard deviation A(y). In addition, we also consider a step function for the 
initial condition, i.e., 

c(r, 0) = c 6(L y /2 - y) - c 6(y - L v /2) + £(r), (4) 

where 6 is the step function. Such an initial condition would be appropriate if the system 
is allowed to phase separate and equilibrate and then quenched to another point inside the 
coexistence curve. 

In Refs. [[J and ||, it was assumed that the only effect of the concentration gradient was 
to change the local average concentration, c(y). The fact that the region in which phase 



separation has occurred grows with time was interpreted to be because the timescale for 
phase separation increases with increasing deviation from the critical concentration |J . That 
is, inside the spinodal, the linear growth timescale diverges as c(y) approaches the mean field 
spinodal while, outside the spinodal, the nucleation timescale diverges as c(y) approaches 
the coexistence value. However, there is also a second possibility. Phase separation in 
the region where c(y) ~ (near the critical concentration) may serve as a seed for phase 
separation in the rest of the system. That is, the time at which phase separation occurs 
is not determined by the local c but by the fact that phase separation has occurred in an 
adjoining region. To test for this possibility we used a linear gradient profile (Eq. (||)) with 
the initial randomness restricted to a small strip in the middle of the system. Specifically, 
A(y) = 0.01 for 0A5L y < y < 0.55L y and A(y) = otherwise. Any phase separation outside 
the middle strip is therefore initiated from this strip. 

Figure p] shows the patterns formed after such a quench with cq = 0.8. Immediately 
after the quench a convoluted structure similar to that for a critical quench is formed in 
the central region. There are well-defined boundaries within which phase separation has 
occurred and outside which it has not. The phase separated region grows leaving behind 
a droplet-like structure |T2[]. This zone continues to grow until it reaches where the local 



concentration is of order the mean field spinodal, c ~ c sp = ±l/\^3. For intermediate times, 
a convoluted interface exists in the middle, separating regions of majority plus phase from 
that of majority minus phase. This interface becomes smoother with time but does not 
become planar during our simulations. Even at t = 6400, it has some curvature correlated 
with the remaining droplets. 

We repeated the experiment with initial randomness over the whole system. Very similar 
behavior was found with the phase separated region expanding at about the same rate as 
before. Therefore we conclude that, even for uniform randomness, the time at which phase 
separation occurs is determined by phase separation in a neighboring region rather than 
there being uniform phase separation (based on the local concentration) over the entire 



system ||13|| . Which of these two mechanisms dominate depends on the relative magnitudes 



of the initial concentration gradient, the initial randomness and thermal noise. Elucidating 
this dependence should prove interesting but will not be done here. 

Figure shows the concentration profile averaged in the x direction, c(y, t) = 
L~ 1 Jq x cIx c(r,t). There is no change in the profile at early times. However, with in- 
creasing time, c(y,t) approaches its bulk equilibrium values of ±1 first at the boundaries of 
the system (where c(y,0) is in the nucleation regime). This equilibrium region invades the 
phase separated regions very slowly since the invasion requires the elimination of droplets. 
The Ising model simulations did not observe this slow change in c(y, t) due to the smaller 
systems and times used ||. 

Significantly different behavior was found for the step function initial condition (Eq. 
(f|)). Figure |3] shows the patterns obtained with Co = 0.5 and A(y) = 0.1 for all y. Although 
the central interface is unstable at very early times due to the flux across it (the Mullins- 
Sekerka instability) ||14|| , it rapidly becomes planar since the flux decreases rapidly with time. 
Droplets then emerge everywhere except in a 'depletion' region near the central interface. 
Both the droplet size R(t) and depletion region thickness, Lu(t), grow in time. A plot of 
c(y, t) (Fig. |j) also show the substantial difference between the two initial conditions. Here 
c(y, t) first approaches the equilibrium interface profile in the central region. The region in 
which c(y,t) is ±1 grows around the central interface and is precisely the depletion region 
seen in the patterns. Simulations with cq = 0.65 > c sp and A(y) = 0.01 were also performed. 
In this case, no droplets were formed because fluctuations are not strong enough and one 
only observes the planar central interface. 

Therefore, for Co < c sp , one of the dominant late stage processes is the growth of droplets. 
However, since the final equilibrium state is the two coexisting phases separated by a planar 
interface, therefore there must also be a macroscopic flux transferring mass across the central 
interface. (In fact, the total area occupied by the droplets decreases even though the char- 
acteristic size is increasing.) Therefore a measure of the approach to the final equilibrium is 
the total mass transferred per length, m(t), 



5 



1 / [Ly fLy/2 \ rL x 

m(t) = — / dy- dy) dx [c(r, t) - c(r, 

L x \JLy/2 JO J Jo 



2 dt'j(t'), (5) 
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where j(t) is the average current density across the middle of the system, 



(6) 

y=L y /2 



This quantity can be measured directly. However, it may be experimentally simpler to 
examine a first order transition in which the motion of the domains are determined by the 
rate at which heat is removed. In this case the thermal current is analogous to j(t) discussed 
above. 

Figure [| shows j(t) for step function initial conditions with Co = 0.5 and Co = 0.65. For 
Co = 0.65 we find j(t) ~ t^ 1 ! 2 for all times. At early times, j(t) for Cq = 0.5 is very similar. 
However, at t ~ 100, droplets begin to appear and there is a rapid decrease in the current 
density. There is a crossover to the late time behavior, j(t) ~ t~ 2 / 3 , which seems to begin 
when the local equilibrium is established. For finite L y , there is a further drop in the current 
density when the last bubble disappears. This last time depends on L y , for L y = 64, the 
final regime begins at t w 3000. 

To understand the behavior of j(t), we consider the chemical potential /x which, in 
contrast to c(r, t), is continuous on all lengthscales. If Co > c sp + A no phase separation 
occurs in the bulk phases. The dynamics are effectively one-dimensional (see Fig. |5]). Except 
at the central interface, we can expand \i as \i = 26c + 0(5c 2 ) where c = ±1 + 5c. (The 
sign depends on the which side of the interface one is on.) The change in c(y, t) occurs first 
near y = so the appropriate boundary conditions are fi(y = ±oo, t) = i/i^ = T c o(l ~ c o) 
while, local equilibrium gives = 0,t) = 0. Substituting 5c ~ /i/2 into the Cahn-Hilliard 
equation gives, 

dt dy 2 ' { ] 



Therefore this regime is diffusion controlled [|l^]. The asymptotic current density, j(t) 



—dn(y,t)/dy\ y= o, is then j(t) = —[too/ V^t in agreement with our numerical results (for c 
near 1). 

The above analysis no longer holds if patterns are present. The dynamics are no longer 
one-dimensional and the droplets are important (see Fig. |5|). Figure |] shows that most 
of the variation in c(y,t) occurs within the depletion width, Lo(t), while in the bulk, the 
patterns are independent of y. The growth of the droplets in the bulk phases is that of 
an off-critical quench with R ~ i 1 / 3 |4],[5|, . Therefore, jx^ must be time dependent 



with local equilibrium giving fj,^ ~ 1/R ~ ^t" 1 / 3 [|l(J. To complete the picture we need 
the behavior of Lr>(t). In principle, Lr>{t) may grow with a dynamic exponent different 
from R(t). However, results for the Cahn-Hilliard equation near a wall shows the exponents 



are the same in both parallel and perpendicular directions [1/]. Therefore we assume that 



-^d(^) ~ R{t) ~ t 1 ^ 3 and the asymptotic current density becomes 

*> ~ A ~ (8) 

in agreement with our numerical results. 

If our assumptions hold, the chemical potential averaged over x has the scaling form 

f i{v,t)=t- 1 "g{v/t 1 t*). (9) 

Figure ^| shows the scaling function, g(z), for times from t = 400 to t — 6400. In agreement 
with our discussion, most of the variation in fi(y, t) is in the depletion region and approaches 
time-dependent bulk values outside it. The collapse of g(z) is reasonable indicating that the 
scaling form of fJ,(y,t) is correct and that our arguments for j(t) ~ t~ 2 / 3 are valid. 

In summary, we examine phase ordering dynamics for constant gradient and step func- 
tion initial conditions. We demonstrate an additional mechanism for the expanding phase 
separated region observed in experiments. We discuss how the concentration profile changes 
with time. We measure the flux characterizing the evolution to the final equilibrium state. 
Its asymptotic behavior is obtained and understood in terms of scaling arguments. 
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FIGURES 

FIG. 1. Concentration patterns for a constant gradient initial condition with Co = 0.8, 
L x = 256, L y = 512 and A(y) = 0.01 for 0A5L y < y < 0.55L y . 

FIG. 2. The concentration profiles for parameters in Fig. |l] (averaged over five configurations). 
The equilibrium values c(y,t) = ±1 are first reached at the boundaries. The oscillations at late 
times are because the bubbles arrange in rows next to the line where c(?/, 0) — c S p. 

FIG. 3. Concentration patterns for a step function initial condition with cq = 0.5, A(y) = 0.1 
and L x = L y = 512. 

FIG. 4. The concentration profiles for the parameters in Fig. |3| (averaged over 5 configura- 
tions). 

FIG. 5. j(t) vs. t in one and two dimensions for step function initial conditions. The solid line 
is j(t) ~ r 2 / 3 . For c = 0.65, j(t) ~ t" 1 / 2 . 

FIG. 6. Scaled chemical potential, g(z) = t 1 ^ 3 ^(zt 1 ^ ,t), for the parameters in Fig. ^. The 
unsealed /j, is shown in the inset. 
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